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1. Introduction 

GALPROP is a software package for numerical Galactic cosmic-ray propagation and related 
emission processes. Its origins go back more than 20 years (1; 2; 3). The first mature version of 
GALPROP was presented in (4; 5). and it has been developed further in the following years. It is 
described in a wider context in (6). A recent related review is (7). 

The extensions described here were made based on the C-i-i- version (v54) publicly released 
four years ago in 2011*. These developments were prompted by new requirements and features 
which appeared desirable, including suggestions by users. The extended version is available^; 
it has been downloaded over 400 times since its initial release in 2013, showing the community 
interest. It is an ongoing development, with regular updates. GALPROP is maintained under the 
configuration control system subversion (svn). I briefly describe each of the new features, grouped 
according to topic, with desiderata for future developments. Details can be found in the Explanatory 
Supplement document included in the package. 

2. Physical processes 

2.1 Hadronic production cross-sections 

There has been progress in hadronic y-ray production physics (8; 9; 10) and this is reflected in 
new options including combinations appropriate to low and high energy protons^. A correction for 
nuclei heavier than Helium in both CR and interstellar gas is provided. Corresponding options for 
secondary positrons and electrons are foreseen in future. 

2.2 Hadronic energy losses 

Energy loss of nuclei by hadronic processes (pion production) was not included previously, 
only ionization and Coulomb losses were present, so that energy conservation was not respected 
since there is hadronic production of y-rays and leptons. Hadronic losses are comparable to ioniza¬ 
tion and Coulomb losses around 1 GeV. This is a potentially important process especially in regions 
of high gas density. It has been included in a simple treatment based on (11; 12) for protons and 
Helium. This is a complex problem however, and further developments are desirable. In particular 
how to treat losses in the context of nuclear fragmentation, beyond the approximation of conserved 
energy per nucleon, is difficult. 

2.3 Deuterium production by pp fusion 

At energies up to 1 GeV, deuterium production by pp fusion is significant and has been in¬ 
cluded^ using the formulation in (13). See also Picot-Clemente et al. this conference, paper 
PoS(ICRC2015)555. 

'http://galprop.stanford.edu. Last update 2011. 

^https://sourceforge.net/projects/galprop 

^Thanks to Chuck Dermer, Michael Kachelriess and Sergey Ostapchenko for providing the relevant code, and much 
help on this topic. 

^Thanks to Nicolas Picot-Clemente for suggesting this topic, and help in its implementation. 
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2.4 Synchrotron radiation 

The treatment of synchrotron radiation from electrons and positrons has been updated to in¬ 
clude polarization using the regular component of the magnetic field. For a detailed description 
and results see (14; 15). GALPROP models of the synchrotron spectral index have been used in 
the analysis of Galactic microwave emission observed by Planck (16). 

2.5 Free-free emission and absorption 

Although not a CR-related process, radio surveys include free-free emission from thermal 
electrons, so skymaps are now produced. These are relevant at frequencies above 1 GHz, and 
dominated for WMAP and Planck. The model is based on the NE2001 thermal electron model. 
Free-free absorption is important at radio frequencies below 100 MHz, and this can be applied to 
both synchrotron and free-free emission. Again NE2001 is used. See (15) for details and results. 
Developments beyond NE2001 would be desirable in future since this has known problems. 

3. Cosmic-ray propagation 

3.1 Accurate solution options 

The original numerical scheme for solving the CR propagation equations uses a Crank-Nicolson 
scheme with operator splitting, and a procedure to accelerate the solution by varying the time step 
from large to small values. Both the operator splitting and time-reduction schemes have the ef¬ 
fect that the true solution to the steady-state equation is not obtained (see Explanatory Supplement 
and also (17)). This can be checked by detailed diagnostics of the solution, provided as a GAE- 
PROP option. In order to obtain accurate solutions, an explicit scheme is now available, which 
uses small enough time-steps that none of the mentioned approximations are required, and demon¬ 
strably converges to the steady-state solution. Detailed analysis is in the Explanatory Supplement. 
The disadvantage of this method is the long computational time on account of the small time steps. 
This is clearly a brute-force approach, and the application of modern numerical methods, such as in 
PICARD (18; 17), which can solve the equations fast and accurately, would be preferable in future. 

3.2 Anisotropic diffusion 

Spatial diffusion was assumed to be isotropic, but theory predicts that it is anisotropic, being 
slower perpendicular to the regular field, hence slower in the z-direction. An option for separate 
diffusion coefficients in the Galactic plane and perpendicular to it has therefore been implemented. 
In future, a more specific relation to the regular B-field would be desirable. 

3.3 Convective transport 

Convection (Galactic wind) was originally represented as a velocity increasing linearly from 
the Galactic plane, which is rather unphysical. A more plausible form has been introduced, which 
still respects the condition of zero velocity at z = 0 (since the sign must reverse there) and which 
smoothly increases to a constant value beyond some z. The form (using a tank function) is con¬ 
trolled by two parameters for the strength and z-dependence of the wind. In future a physical 
description of the wind would be desirable. 
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3.4 Boundary conditions 

The spatial boundary condition imposed has been zero CR density at all boundaries in 2D 
(/?,z) or 3D {x,y,z)- This is clearly not physical since even with free-escape, the density is non-zero 
there. The condition was imposed by simply forcing a zero value after each time-step, rather than 
acually solving with this condition. As an alternative, the boundary condition is not imposed, so 
particles flow out but not back in at the boundary. The density at the boundary is non-zero. A more 
physical approach actually handling the physics of the boundary, would be desirable in future. 

3.5 Primary positrons 

There is good experimental evidence for a primary component of positrons (PAMELA, AMS- 
02), in addition to secondaries from hadronic interactions. Previously only secondary positrons 
were included. Primary positrons can now be included, with an injection spectrum and spatial 
source distribution independent of primary electrons. Note that primary positrons also contribute 
to /-rays via inverse-Compton and bremsstrahlung, and to synchrotron emission, so are important 
to take into account. 

3.6 Injection spectrum 

Additional spectral breaks in the injection spectra of nuclei, electrons and positrons are now 
available, for more flexibility in fitting observations. 

4. Galactic structure 

4.1 Magnetic fields 

Additional models of the Galactic magnetic field, including the regular field described by 
(19; 20), are now included. 

5. Format and examples 

5.1 HEALPix 

Skymap output (/-rays , radio) is provided in HEALPix (21) as the format of choice for full- 
sky coverage with uniform pixel size. New versions of HEALPix have become available^, and the 
GALPROP code has been adapted to enable their use. The format has been adjusted so that the 
energy (or frequency for radio) is stored as PITS columns rather than the original format of vectors 
within pixels. This is compatible with the Aladin visualization package®, and also more convenient 
for reading by user software. 

5.2 Examples and reference output 

A variety of example parameter files (galdef files) are provided to illustrate the many features 
including the new ones. As an installation check, a full set of output files is also provided for one 
sample case. 


^http://healpix.jpl.nasa.gov 

®http://aladin.u-strasbg.fr 
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5.3 Bugfixes 

Various (minor) bug fixes to the original version were made during the development. 

6. Related software 

A plotting package GALPLOT, for cosmic rays, /-rays and radio, compatible with GALPROP 
output, is now available^. It is compatible with the cosmic-ray database CRDB (22), and can handle 
Fermi-LAT and other /-ray data and various synchrotron surveys. GALPLOT can also be used for 
Galactic source population synthesis, as described in section 6 of (23). 

Separate routines for computing synchrotron^, anisotropic inverse Compton scattering, and 
models for inverse-Compton emission from solar and stellar heliospheres are also available^. 

7. Outlook 

Other projects complementary to GALPROP are in progress. DRAGON (24) includes various 
new features. PICARD (18; 17) allows full 3D models with high resolution using advanced nu¬ 
merical methods. The semi-analytical USINE (25) package is becoming public (see Maurin, this 
conference, paper ID 296). These activities will stimulate further GALPROP developments. 
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